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ABSTRACT 

We use a sample of 262 spectroscopically confirmed star-forming galaxies at redshifts 2.08 < z < 2.51 
to compare Ha, UV, and IR star-formation-rate diagnostics and to investigate the dust properties of 
the galaxies. At these redshifts, the Ha line shifts to the AT s -band. By comparing /\ s -band photometry 
to underlying stellar population model fits to other UV, optical, and near-infrared data, we infer the Ha 
flux for each galaxy. We obtain the best agreement between Ha- and UV-based SFRs if we assume 
that the ionized gas and stellar continuum are reddened by the same value and that the Calzetti 
attenuation curve is applied to both. Aided with MIPS 24 /jm data, we find that an attenuation 
curve steeper than the Calzetti curve is needed to reproduce the observed IR/UV ratios of galaxies 
younger than 100 Myr. Furthermore, using the bolometric star-formation rate inferred from the UV 
and mid-IR data (SFRir+SFRuv), we calculated the conversion betwee n the Ha lu minosity and SFR 
to be (7.5 ±1.3) x 10 -42 for a Salpeter IMF, which is consistent with the iKennicuttl (Il998l ) conversion. 

The derived conversion factor is independent of any assumption of the dust correction and is robust 
to stellar population model uncertainties. 

Subject headings: galaxies: evolution — galaxies: higli-redshift - galaxies: star formation 


1. INTRODUCTION 

One of the most important diagnostics in understand¬ 
ing the evolution of galaxies is the star-formation rate 
(SFR). The evolution of the SFR of galaxies can give 
clues as to how galaxies were enriched with heavy el¬ 
ements, how they build up their stellar mass through 
cosmic time, and helps us to understand the bolomet¬ 
ric output of galaxies. At redshift z ~ 2, when the 
universe was just ~ 3 Gyr old, star-formation activity 
in the universe was at its peak and galaxies were in 


the process of assemblin 

K most of their stellar mass 

(see Reddv & Steidel!2009 

:lBouwens et al.ll2010; ShaDlevI 

2011; IMadau & Dickinson 

2011. Studying this critical 


epoch is essential to gaining a better understanding of 
the evolution of the progenitors of the local galaxy pop¬ 
ulation. 

The ultra-violet (UV) continuum (1500 to 2800A) in¬ 
tensity of a galaxy is one of the most commonly used 
diagnostics for the SFR as it is observable over a wide 
range of redshifts and intrinsic luminosities. It is sensi¬ 
tive to massive stars (M* > 5 M 0 ), making it a direct 
tracer of current SFR. By extrapolating the formation 
rate of massive stars to lower masses, for an assumed 
form of the init ial mass funct ion (IMF), one can estimate 
the total SFR (IMadau et al.lll998h . Another widely used 
diagnostic for measuring the SFR is nebular emission, 
with Ha being the most common because of its higher 
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intensity compared to the other hydrogen recombination 
lines such as H/3, Paa, Pa/3, etc., and it is easier to inter¬ 
pret than the Lya line. Ha is an “instantaneous” tracer 
of SFR because it is sensitive only to the most massive 
stars (M* > 10 M 0 ). However, it becomes more chal¬ 
lenging to observe Ha from the ground at z > 1 because 
the line is redshifted to the near-IR where the terrestrial 
background is much higher than at optical wavelengths. 

The main disadvantage of using UV/optical luminosi¬ 
ties as tracers of the SFR is their sensitivity to dust at¬ 
tenuation. The dust absorption cross-section is larger for 
shorter wavelengths and choosing the appropriate atten¬ 
uation curve to correct the observed luminosities plays an 
important role in determining intrinsic physical quanti¬ 
ties. Aside from the assumed attenuation curve, the ge¬ 
ometry of dust with respect to the stars can lead to differ¬ 
ent color excesses, E(B—V ), between the ionized gas and 
the stellar continuum. E{B — V ) is the color excess mea¬ 
sured between the B and Vbands, E(B — V) = As — Ay, 
where A\ is the total extinction at wavelength A in mag¬ 
nitudes. In particular, the nebular recombination lines 
arise from the HII regions around the most massive O 
and early-type B stars (with masses of M* > 10 M 0 
and main sequence lifetimes of < 10Myr). On the other 
hand, for a Salpeter IMF, solar metallicity, and a con¬ 
stant or rising star-formation history, the UV continuum 
in starburst galaxies originates from stars over a broader 
range of mass tha t inclu d es later -type B stars w ith life¬ 
times < 100 Myr (IKennicuttl 1 19981 IMadau fc Dickinsonl 
I2014T) . These older non-ionizing stars have more time to 
migrate to regions of lower dust density in the galaxy, 
while H-ionizing stars with shorter lifetimes do not have 
enough time to escape from their dusty birthplace or let 
the parent molecular clouds to dissipate. As a result, 
the nebular lines can be subject to a higher degree of 
r eddening than t he UV continuum. 

iCalzetti et al.l (119941) found that the nebular emission 
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is more attenuated than the stellar continuum at the 
same wavelengths fo r a sam ple of local UV-bright galax¬ 
ies. Subsequently. iCalzetti et all (|2000D studied a similar 
sample of local galaxies and argued that the color-excess 
is 2.27 times larger for the nebular emission lines than for 
the stellar continuum. This relationship was derived un¬ 
der the assumption that the Calzetti curve is applied to 
the stellar continuum and a Galactic extinction curve is 
app lied to the nebular emission lines. In a separate study, 

ICid Fernandes et~ahl ()2005l ) found that for an SDSS sam¬ 
ple of 50000 galaxies, .4 v n einti a-r ^ 1.8 Ay stellar i both 
assuming a lCardelli et all ( 1989 1 extinction curve. 

There have been several studies of relative attenuation 
of the stellar continuum and nebular regions in high- 
redshift galaxies. In most of the high-redshift studies, 
the Calzetti curve is assumed fo r both the nebular and 
stellar regions. Frequently, the ICalzetti et al. (1200 0) re¬ 
lation between the stellar and nebular E(B — V) is used 
to correct the nebular lines, while this relation was origi¬ 
nally derived based on using a Galactic extinction curve 
for the nebular lines. At present it is unclear what type 
of attenuation curve should be applied to the nebula r re- 
gions. These issues are discussed in detail in ISteidel et all 
(120141) . 

At z ~ 2, IForster Schreiber et all (I2009T) selected a 
sample from various imaging surveys in the optical, 
near-IR, mid-IR, and sub-millimeter regime, with SED- 
derived SFRs of ~ 10 - 1000 M 0 yr -1 and stellar masses 
> 10 10 M 0 , assuming a Chabrier IMF. In their study, the 
predicted Ha luminosities from the best-fit SED mod¬ 
els for 62 galaxies were compared with the observed line 
luminosities measured from integral field spectroscopy. 
The latter were then corrected for reddening assuming 
the Ay return ed from stel l ar popula tion modeling of 
their galaxies. IForster Schreiber et all ( 2009 ) found that 
assuming a factor of two redder color excess towards the 
HII regions relative to the stellar continuum yielded the 
best agreement between the predicte d and observed Ha 
luminosities. Similarly, IWuvts et all (1 201 ll) investigated 
25 star forming galaxies at z ~ 2 and compared the Ha 
SFRs, dust-corrected by the Calzetti attenuation curve, 
with the SED modeled SFRs. They showed that Ha 
SFRs are consistent with the SED SFRs, provided that 
E(B - V) neb = 2.27 E(B - V) ste iiar when assuming the 
C alzetti cu rve fo r bo th. 

iMuzzin et al.l (f2010h presented two galaxies at red- 
shifts z = 2.122 and 2.024, and added ten othe r galax ¬ 
ies at similar reds hifts f rom Ivan Dokku m et all (|2005l ). 
iKriek et~ahl (I2007T) , and IKriek et all (120091 ) to compare 

A-y,nebular derived from the Balmer decrement () 
with A y,stellar derived from the SED models, both as¬ 
suming the Calzetti curve. The only two galaxies with 
detected H/3 line measurements showed evidence for 
higher nebular attenuation, but about half of the non¬ 
detections were consistent with equal attenuation for th e 
nebular and stellar regions. In lYoshikawa et ahl ( 2010 ). 
Ha attenuation was investigated through several com¬ 
parisons for a sample of optically-selected star-forming 
galaxies at z ~ 2. Aided with n e ar-inf rared spectroscopic 
observations, lYoshikawa et al.1 (120101) used the Balmer 
decrement along with the Calzetti curve to correct the 
Ha SFR. The dust-corrected SFR was then compared 
with the Ha SFR corrected by the SED color-excess and 


also with the SFRs inferred from UV, IR, radio, and X- 
ray. The comparisons yield contradictory results. Their 
analysis favored a higher nebular color-excess for galaxies 
with larger intrinsic S FRs. 

In s eparate studies. lErb et, al.l (|2006d ) and lReddv et all 
(2010) used a large sample of UV-selected galaxies at 
z ~ 2 and corrected both the Ha and UV SFRs with 
the color-excess derived from the best-fit SED models 
and applied the Calzetti attenuation curve. They showed 
that using different color excesses for the continuum and 
lines generally resulted in Ha-inferred SFRs that over¬ 
predicted those calculated from the dust-corrected UV 
emission and direct IR measurements. 

The apparently conflicting results mentioned above 
may be reconciled if the relation between the nebular 
and stellar color-excesses depends on t he physical prop¬ 
erties of galaxies, such as their SFRs (lY oshikaw a et all 
20101) or specific SFRs ()Wild et al.l 1201 l b For example, 
Wild et all (1201 ll) showed that the ratio of the line opti¬ 
cal depth to the continuum optical depth decreases with 
increasing sSFR and at sSFR~ 10~ 9 yr ~ 1 the li ne-to- 
contin uum optical depth ra tio reaches t h e ICalzetti et all 
(|2000j) ratio assuming the iWild et al.l (|2011l ) attenua- 
tion curve. iPrice et all (|2014D investigated the atten¬ 
uation of the nebular regions at z ~ 1.5 using the 
Balmer decrement from stacked HST grisrn spectra as¬ 
suming the Calzetti curve, and found no stron g trend of 
A v , neb/ A v , stellar with SFR or sSFR. In lPrice et all (12014D 
study, at sSFR~ 10~ 85 yr _1 the ratio is consistent with 
1 . 

The primary goal of this paper is to understand the 
relationship between the UV and Ha emission in high- 
redshift star-forming galaxies, with a large sample that 
is immune to uncertainties in slit los s correct ions that 
affects the Ha flux estimation (c.f. lErb et all l2006d : 
lYoshikawa et al.ll2010 ), the small samp le sizes inherent in 
previous spectrosco pic studies of Ha (|Kriek et al.l 120071 : 
IMuzzin et al.l l2010h . and not subject to the bias of se¬ 
lecting high equival ent w idth objec ts from narrow-band 
selected samples dGarn et al.ll2010D . We consider in our 
analysis only spectroscopically-confirmed galaxies, en¬ 
abling us to estimate the Ha flux from broad-band pho¬ 
tometric excess techniques without the additional uncer¬ 
tainties that plague photometric redshifts (IWuvts et al.1 

mi ii). 

The impact of nebular lines on the broad-band 

photometry was known an d stu died_for many 

years (e.g.. iGuiderdoni fc Rocca-Volmerangel 1987 : 
Fioc fc Rocca-Volmerangel 119971 : iZackrisson et al.1 [2008; 
de Barros et ah 120141) . Using the photometric excesses 
to determi ne the line strengths was pioneered by 
IShim et all ( 2011 ). where they showed that the excess 
in Spitzer/TRAC 3.6 /jm relative to the SED model 
continuum is due to the redshifted Ha emission line for 
t heir sample of galaxies at 3.8 < 2 < 5.0. iStark et ahl 
(2013) also investi gated a sample of g alaxies at the same 
redshift range of iShim et al.l (120111) and inferred the 
Ha emission line strengths by comparing the observed 
flux in Spitzer /IRAC 3.6 pm band and the continuum 
flux as expecte d jrom the SED model. Following that, 
iSchenker et al.l (|2013f ) verified the photometric excess 
technique by applying it to a small sample of 9 galaxies 
at 3.0 < 2 < 3.8; for 8 galaxies the [O hi] line fluxes 
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inferred by the same technique as iStark et ~ahl (|2013T ) 
were within a factor of 2.5 of the spectroscopically 
measured [O hi] fluxes. 

An additional advantage of our study is that we include 
IR data to independently assess the dust-obscured SFR. 
Comparing Ha, UV, and IR-inferred SFRs in a statistical 
sense allows us to understand how to correct extinction- 
sensitive measures of SFR for the effects of dust. 

The outline of this paper is as follows. In j[2]we discuss 
the properties of our sample, the assumptions that have 
been made to model the stellar populations using the 
rest-frame UV to near-IR photometry, the Spitzer /MIPS 
photometry and stacking method. A detailed descrip¬ 
tion of how we estimated the Ha line flux is provided 
in El The analysis of the MIPS 24 pm data and IR lu¬ 
minosities is described in El In El we compare the two 
different tracers of SFR - Ha and UV - and discuss the 
dust correction recipe most consistent with our measured 
values. E] focuses on combining SFR diagnostics (e.g., 
Ha and UV with IR) to deduce bolometric SFRs. The 
results a r e sum marized in El Throughout this paper, a 
iSalpetef ( 195 5) IMF is assumed and a cosmology with 
H 0 = 70 km s -1 Mpc -1 , Ha = 0.7, and H m = 0.3 is 
ado pted. All magnitudes are given in the AB system 
(iOke fc Gunnl[l983 ). 

2. SAMPLE 

2.1. Sample Selection and Optical Photometry 

The galaxies used in this study are drawn from a 
parent sample that is part of an imaging and spec¬ 
troscopic survey o f UV-selected galax ies at z ~ 2-3 
(ISteidel et all 120041 : iReddv et al.ll2012bll . The galaxies 
were selected based on the B X, BM, and Lyman break 
galaxy rest-UV color c riteria (jSteidel et all 120031.12004 
lAdelberger et all 12004 1. where U n GR optical data were 
obtained with the Palomar Large Format Camera (LFC) 
or Keck Low Resolution Imaging Spectrograph (LRIS; 
iSteidel et al.ll2003l 120041 . Rest-UV spectroscopic follow¬ 
up with Keck/L RIS was conducted f or gal axies brighter 
than R = 25.5 (ISteidel et al.ll2003l l200l . Near-IR J 
and K s imaging was obtained using the Palo mar /WIRC 
and Magellan/PAN IC instruments dShaplev et al .1 [20051 : 
IReddv et al.1 l2012bl) . A-band (F160w) data were ob- 
tained with the Hubble Space Teles cope WFC3 camera 
(|Law et alJl2012t IReddv et al.ll2012bh . All galaxies in the 
sample have coverage in at least one of the Spitzer/IRAQ 
four channels (3.6, 4.5, 5.8, and 8.0 ^m; R.edd v et al.1 
l2006al l2012bl) . The objects are located in the GOODS- 
North field and 11 additional fields that were primarily 
selected to have one or more relatively bright background 
QSOs for studying the interface between the i ntergalac- 
tic medium (IGM) and galaxies at z ~ 2-3 (ISteidel et al.1 

[2001 [M3). 

Out of the final sample of 2283 objects with spectro¬ 
scopically confirmed redshifts, a subset of galaxies is se¬ 
lected based on the following criteria: (1) the galaxy 
must be covered by the A-band imaging, (2) it must have 
a redshift in the range 2.08 < z < 2.51 so that the Ha 
line falls into the K s band, and (3) it must be detected 
in at least two of the IRAC channels or one of the IRAC 
channels and either J or F160w bands. The third condi¬ 
tion ensures a more robust estimate of the stellar mass 
and the continuum level at 6564A. Furthermore, AGNs 


(making up « 9% of the parent sample) were identified 
by either strong UV emission lines (e.g., Lya, CIV) or by 
a power law SED through the IRAC bands. These AGNs 
are removed from our sample. Eventually, 262 galaxies 
remain that satisfy the aforementioned criteria. 

2.2. Stellar Population Modeling 

For each galaxy in our sample, the best-fit stellar popu¬ 
lation model is found by using the rest-frame UV through 
near-IR broad-band photometry. As mentioned above, 
all the galaxies in our sample have spectroscopically- 
confirmed redshifts, thus removing a key degeneracy in 
the modeling of the stellar populations. In addition, for 
a better estimation of stellar mass and age, all galaxies 
in our sample have at least two detections long-ward of 
the 4000 A break - excluding A' s -band. 

C. C hariot & G. Bruzual (2007) models with a lSalpeteil 
(11955T) IMF and solar metallicities are used for the fit¬ 
ting. For each individual galaxy, different star-formation 
histories are assumed, including constant, exponentially 
declining, and exponentially rising with characteristic 
timescales of r = 10, 20, 50, 100, 200, 500, 1000, 2000, 
5000 Myr for exponentially declining and r = 100, 200, 
500, 1000, 2000, 5000 Myr for exponentially rising his¬ 
tories. Ages are allowed to vary between 50 Myr and 
the age of the universe at the redshift of each galaxy. 
The > 50 Myr limit corresponds to the typical dynam¬ 
ical timescale of star-forming galaxies at z ~ 2 as in¬ 
ferred from velocit y dispersio n an d size measurements 
of these galaxies (IReddv e t al.ll201 2bjl. Fo r interstellar 
dust obscuration, the ICalzetti et al.[(j2000t l attenuation 
curve is used, allowing E(B — V) to vary between 0.0 
and 0.6. The x 2 values have been determined for each 
set of observed broad-band and model magnitudes. The 
best-fit model is determined through y 2 minimization. 
There is generally no significant difference between the 
best-fit x 2 values of the six different population models 
(constant, exponentially rising, and exponentially declin¬ 
ing star-formation histories, for each considering all ages 
and ages greater than 50 Myr). As previous studies have 
shown, the assumption of declining star-formation his¬ 
tories at these redshifts results in systematically lower 
SED-i nferred SFRs compared to the observe d IR+UV 
SFRs dWuvts et all 20 lit Reddv et alJl2012bll . Further¬ 
more, IReddv et all (12012bh showed that assuming a con¬ 
stant star-formation history for z ~ 2 galaxies predicts 
specific SFRs (SFR/M*) at higher redshifts (z ~ 2.6) 
that are substantially larger than the observed values. 
Given these, we adopt the models that assume expo¬ 
nentially rising star-formation histories with ages greater 
than 50 Myr. 

The SED models used to fit the observed magnitudes 
do not include nebular emission lines. For example, the 
Ha line can significantly affect the photometry and we 
use its contribution to the A-band to estimate the Ha 
line flux. The [O hi] emission line is the other strong line 
which falls into the F160w filter given the redshift range 
of our galaxies. Only 20% of the galaxies have F160w 
observations, for which we did not correct the broadband 
photometry for the contamination. The SED-inferred 
SFR of these galaxies is consistent with the SFR(UV) 
estimates within the uncertainties. At this redshift, J 
band is contaminated by the [O ii] emission line, but this 
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line is generally weaker than the Ha line and its effect 
on the SED inferred SFRs is negligible compared to the 
uncertainties^] 

2.3. MIPS data 

To further investigate the bolometric properties of our 
sample, we use Spitzer /MIPS 24 rtm wherever available . 
Out of 12 fields, GOODS-North (iDickinson et, al.l 12003(1 . 
and four other fields (Q1549, Q1623, Q1700, and Q2343; 
Reddy & Steidel 2009) have MIPS 24 (im coverage to a 
typical 3 ct depth of 10-15/iJy. 

Photometry on 24 /jm images is performed by using 
point-spread function (PSF) fitting with priors deter¬ 
mined by the locations of the objects in the higher res¬ 
olution IRAC images (IRAC data exist in all fields). A 
40x40 pixel region centered on each target is extracted 
with pixel size of l'/2. PSFs are then fitted simultane¬ 
ously to all known sources in the sub-image and one ran¬ 
dom background position. This procedure is repeated 
many times to obtain sufficient statistics for proper back¬ 
ground estimation based on the random background flux 
measurements. The other source of uncertainty is Pois¬ 
son noise, which for objects in our sample is negligi¬ 
ble compared to the background dispersion. We remove 
objects whose photometry may be compromised due to 
blending with nearby sources, using the criteria specified 
in lReddv et al.l (j2010H . This results in 115 galaxies with 
secure PSF fits. 

Out of 115 objects with MIPS data, 47 are detected 
with S/N <3. Undetected objects, those with S/N <3, 
are considered using either survival analyses or through 
stacking of the 24 /an data. 

2.3.1. Stacking Method 

Throughout the paper, we employ a stacking method 
to determine the median 24 /mi fluxes of objects that are 
individually detecte d and undete c ted, followin g the pro¬ 
cedures described iu lReddv et al.l (t2010ll2012all . We per¬ 
formed aperture photometry on the stacked images, and 
applied an aperture correction based on the 24 /tm PSF. 
The average background level and noise were determined 
by placing many apertures of the same size used for the 
stacked signal at random positions in the stacked image 
and measuring the average flux level and dispersion in 
flux of these “background” apertures. Furthermore, we 
used bootstrap resampling simulations to estimate the 
dispersion in the fluxes of objects contributing to each 
stack. This was accomplished by creating 100 samples of 
random images in each bin and measuring the standard 
deviation of the median stacked fluxes. The intrinsic dis¬ 
persion in the stacked flux is larger than the background 
error by a factor of ~ 2. 


G00DS-BX1330 



Q1623-BX293 



Mum) 


Q1700-BX575 



X(,u.m) 

Figure 1. Best-fit SED models for three objects in our sample: 
BX1330 in GOODS-North (top), BX293 in Q1623 (middle), and 
BX575 in Q1700 (bottom). All the models assume an exponentially 
rising star-formation history, with a lower limit in age of 50 Myr. 
The K s band data are shown with an open circle in each plot. The 
green and red lines show the best-fit models when we exclude and 
include, respectively, the Tf-band photometry in the fitting. The 
two fits are similar because the K s photometry is down-weighted in 
the fitting due to larger measurement errors. The vertical orange 
line shows the position of the redshifted Ho line. The dashed area 
represents the K s filter coverage. At redshifts of z ~ 2.08 — 2.50, 
the K s filter is affected by the Ho emission line, and the difference 
between the A s -band flux and the underlying stellar population 
model can be used to estimate the line flux. 

by the K s filter transmission curve. The difference be¬ 
tween the observed iv s -band magnitude and the contin¬ 
uum magnitude is used to extract the Ha line flux as 
follows. The observed K s magnitude is considered as 
the sum of the continuum and the Ha fluxes, while the 
SED-inferred K s magnitude represents only the contin¬ 
uum flux. Assuming a Gaussian form for the redshifted 
Ha line, the change between these fluxes will yield the 
Ha flux (Figure [TJ: 


3. Ho AND UV LUMINOSITY 
3.1. Ha Fluxes and Equivalent Widths 

Using the procedure discussed in N2.2I we fit all the 
photometry, excluding AVband, in order to determine 
the continuum level at 6564 A. The continuum magni¬ 
tude is calculated by multiplying the best-fit SED model 

6 For galaxies with similar SFRs at z ~ 2, the typical [O II] line 
flux is ~ 6 X 10“ 17 erg s -1 cm -2 <Kriek et al. ||201~4T ). The ratio of 
the [O II] flux to the typical J-band flux errors in our observations 
is only ~ 0.07. 


If /HaW T K (V) 
I^T k (v) 


(^obs+48'6) (^SEP+48.6) 

10 - 2 - 5 — 10 — 2.5 


( 1 ) 


Here, fHaip) is the flux density of the Ha line in units 
of ergs -1 cm -2 Hz -1 , Xk(f) is the K s filter transmis¬ 
sion curve, and m 0 b s and msED are the observed and 
the continuum magnitudes, respectively. The Ha flux is 
corrected for contamination from the [Nil] lin e based o n 
the mass-dependent [Nn]-to-Ha flux ratios of lErb et al.1 
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Table 1 

Ha Flux Corrections for [Nil] 
Contamination 51 


Mass Range (10 10 Mg) b 

N2 C 

< 0.88 

< -1.22 

0.88 - 2.00 

—i no - * -0 07 

l.UU—009 

2.00 - 3.69 

-0.85l°;°5 

3.69 - 6.03 

—0.78t°;°s 

6.03 - 8.82 

U.OO —o 04 

> 8.82 

U.DD_ 0 02 


a Based on lErb et all (l2006aD 
b Assuming a lSalpeterl d 19551) IMF 

C N2 = log(^ji) 
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(l2006aT) . The stellar masses are determined from the SED 
models and the corresponding [Nil] line contamination, 
as listed in Table [lj is used to corr ect the Ha flux. The 
SFRs are then calculated using the IKennicuttl (U998f) re¬ 
lation to convert the Ha line luminosity to an SFR. 

There are three main sources of uncertainty in the de¬ 
rived Ha fluxes. The largest uncertainty is the pho¬ 
tometric error. The typical error in the observed K s - 
band magnitude is cs 0.27. The second source of un¬ 
certainty is the_[Nii] correction. The uncertainty on the 
lErb et all d2006al) [Nu]-to-Ha line ratios (see Table [1]) is 
added in quadrature to the photometric error. The third 
source of uncertainty is the error associated with the con¬ 
tinuum flux. In order to account for this uncertainty, 
we estimated the continuum flux at 6564A from the 
best-fit model assuming six different star-formation his¬ 
tories: exponentially rising, exponentially declining, and 
constant, for each considering all ages and ages greater 
than 50Myr. For 94% of the galaxies the difference 
in the mean of the continuum fluxes assuming different 
star-formation histories to those assuming a rising star- 
formation history with ages > 50Myr is less than 0.1 
magnitude. The average error in the estimated contin¬ 
uum magnitude is ~ 0.05 mag, which is negligible com¬ 
pared to the observed X s -band magnitude errors. Com¬ 
bining the three sources of uncertainty discussed above 
yields a typical relative error in Ha flux of ~ 0.29. 

The equivalent width of the Ha line is estimated by 
dividing the Ha flux derived from Equation |T] by the 
continuum flux density at the wavelength of the Ha line: 

(l + z)EW 0 = ^. ( 2 ) 

J A 

The fx° nt is the continuum flux density in units of 
erg s -1 cm -2 A -1 that is estimated through the best-fit 
SED model, Fna is the Ha line flux in erg s -1 cm -2 , z 
is the redshift of the galaxy, and EWo is the rest-frame 
equivalent width in A. 

We define whether a galaxy has a “detected” Ha line 
according to the following. Galaxies whose Jf s -band pho¬ 
tometry exceeds the continuum level by more than the 
.Kg-band magnitude error are referred to as “Ha detec¬ 
tions” . Galaxies where the /\ s -band photometry is con¬ 
sistent with the continuum level to within 1 a are re¬ 
ferred to as “Ha non-detections”, and an upper limit of 
1<7 above the measured K photometry is used for these 


Figure 2. Estimated Ha rest-frame equivalent widths as a func¬ 
tion of the error of the observed A s -band magnitude. The circles 
are detected galaxies and are color coded based on their observed 
Ha fluxes (uncorrected for extinction). The red triangles show lcr 
upper limits. Detecting the Ha emission requires that the observed 
K s magnitude lies more that lcr above the underlying continuum, 
where lcr is defined by the K s magnitude error. The K s magni¬ 
tude includes the light both from the Ha emission line and the 
brightness of the continuum at 6564 A. There is a tight correlation 
between the continuum brightness and the K s magnitude error; 
the fainter continuum results in larger uncertainty in the observed 
magnitude. As a result, there are more undetected objects with 
large A s -band errors. On the other hand, the brightest Ha objects 
have bright continuum and small A s -band uncertainties. The de¬ 
ficiency of objects in the upper-left corner is mainly because very 
bright Ha objects that also have bright continuum are rare. 

objects. Out of 262 galaxies in our sample, 149 are de¬ 
tections and 94 are non-detections. There are 19 objects 
with K s magnitudes fainter than the continuum by more 
than lcr. We removed these galaxies from our discussion 
due to their K s photometry being inconsistent with the 
photometry from adjacent bands. 

As discussed in El our method of computing Ha fluxes 
and EWs has the advantage of being immune to slit-loss 
corrections. Our sample covers a wide range of Ha fluxes 
and EWs. Figure [2] shows the distribution of Ha EWs 
and fluxes as a function of the K s magnitude lcr uncer¬ 
tainty. The detection of Ha with the method adopted 
here depends on both the Ha line flux and the bright¬ 
ness of the continuum at 6564 A. As we have defined 
detections to be those objects where the K s magnitude 
exceeds the continuum level by more than lcr , the num¬ 
ber of undetected objects increases for objects that are 
fainter in the continuum (Figure [2]). 

The Ha observed f luxes of 44 objects in common with 
the lErb et al.1 (12006(1) spectroscopic sample are compared 
in Figure [3l The lErb et al.1 (|2006cD spectroscopic fluxes 
are multiplied by a factor of two to account for slit losses 
and the aperture used to extract the spectra (l Erb et all 
I2006d : iReddv et ahll2010H . The slit loss correction factor 
applied to the spectroscopic fluxes depends on various 
factors such as the accuracy of the astrometry, the size 
of the object convolved with the seeing at the time of ob¬ 
servation in comparison with the size of the slit, and the 
accuracy of the mask ali gnmen t during the observation. 
As noted in lErb et all ( 2006 c). the factor of two slit- 
loss correction is an average estimate and the slit losses 
will of course vary from object to object, likely account¬ 
ing for some of the scatter in Figure [3] On the other 
hand, while broadband photometry has its advantages, 
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Ho Flux [erg s" 1 cm' 2 ] _ Erb+06 (multiplied by 2) 


Figure 3. A comparison of observed Ho; fl uxes of t he obj ects in 
common between our sample and that of lErb et al.l ^2006d V The 
diamonds are the detected objects and downward arrows show 1 a 
upper limits, as discussed in the text. The fluxes are measured 
based on the difference between the observed and the continuum K s 
magnitudes predicted from the best-fit model. The error bars are 
dominated by the uncertainty in the observed K s photometry. The 
lErb et ahl H2006cD spectroscopic fluxes are multiplied by a factor of 
2 to account for slit losses. The dashed line denotes the one-to-one 
relationship. 


To estimate the bolometric luminosities of galaxies in 
our sample, we infer the infrared luminosity (integrated 
between 8-1000 pm) as discussed below, and add this to 
the uno bscu red UV luminosity (i.e., observed UV lumi¬ 
nosity; H3.2I) . as in lReddv et ah (1201 Q). 

Infrared luminosities, Lir, are estimated by using the 
Spitzer/MIPS 24 pm observations. The 24 pm band 
is sensitive to the rest-frame 7. 7 pm PAH em i ssion , 
which correlates w i th £tr (e.g.. ICharv fe Elbad 120011 : 
iDale fc Heloul l2002t lElbaz et~ahl 1201 11 1. To convert ob¬ 
served 24/im magnitudes to total IR luminosity (L/r), 
we used several d u st SE P templates, in cl uding those 
of ICharv fc Elbai (1200111 . IDale fe Heloul (1200211 . and 
iRieke et al.1 ( 200T Luminosities determined from the 
rest-frame 8 pm flux density alone, regardless of the dust 
template used, tend to over predict L(IR), particularly 
for LIRGs (e g.. IR eddv et~aLll2012 all and ULIRGs (e.g., 
iNordon et alJ 12010b iMasrnelli et alfboilH . In order to 
account for the luminosity-dependent overestimation of 
the deri ved L(IR1 , we u se a correction equation described 
in lReddv et al.1 (l2012alb The correction is calculated by 
comparing the L(IR) computed from 24, 100, 160 pm, 
and 1.4 GHz fluxes, and the L(IR) determined solely from 
24 pm data for a similarly-selected sample of galaxies: 

logLiR(all data) =0.537 x logLiFt(24/rm) + 5.136. (3) 


our AT s -band measurements suffer from a larger measure¬ 
ment (random) uncertainty relative to the other methods 
given the depth of our ground-based images. Although 
most of the detected Ha fluxes in this study (the black 
symbols in Figure 0|) are higher that the spectroscopic 
fluxes, about half of the sample is estimated as upper lim¬ 
its (the red arrows). Considering both the detections and 
non-detections, the broad-band estimated fluxes are gen¬ 
erally consistent with the spectroscopic fluxes. To quan¬ 
tify the degree of the correlation we conduct a generalized 
Kendall’s tau (r) statistic. The generalized Kendall’s r 
rank correlation coefficient is a non-parametric test that 
permits non-detections in both variables (here, our esti¬ 
mated Ha fluxes). Kendall’s r corr elation coe fficient for 
our flux estimation with the lErb et ahl (120060 1 measure¬ 
ments, excluding the one outlier object, is r = 0.17, with 
significance of p-value = 0.11. 

3.2. UV Luminosities and SFRs 

UV luminosities are calculated using the fluxes of the 
best-fit SED models at 1700 A. The IKennicuttl (I1998T ) 
relation to convert the UV luminosity (over the wave¬ 
length range 1500-2800A) to SFR applies only to galax¬ 
ies where star formation proceeds for > 100 Myr. How¬ 
ever, there are galaxies in our sample with inferred ages 
between 50 and 100 Myr. For these younger galaxies, we 
used an age- and star-formation-history-dependent rela¬ 
tion in the c onversion of the UV luminosity to SFR, as 
discussed in iReddv et al.l (l2012blb The age-dependent 
conversion resu lted i n UV SFRs that are at most 14% 
larger than the IKennicuttl ([1998H SFRs for young galax¬ 
ies. Throughout this analysis, we also use SFRs inferred 
from the best-fit SED models. As expected, within the 
uncertainties, the SED SFRs are highly correlated with 
dust corrected UV SFRs. 

4. THE IR LUMINOSITY AND BOLOMETRIC PROPERTIES 


All luminosities are in Lq. 

Equation[3]is based on the lGharv fc Elba'S : ( 2001 1 mod¬ 
els. There is a factor of ~ 2 variation in L(IR) de¬ 
rived from 24/rm data using IDale fe Heloul (I2002D and 
Rieke et all (12009H dust templates (e.g.. see IReddv et al.1 
20 1 2 all . However, once corrected using the appropriate 
equations similar to Equation^ all the templates result 
in L(IR) estimates that are consistent with each other 
and with the L(IR) computed from 24, 100, 160 pm 1 and 
1.4 GHz fluxes, within the uncertainties of the measured 
IR luminosity as explained in 12.31 (see also IReddv et al.1 
l2012a|h In the subsequent analysis, we use infrar ed lu¬ 
minosities inferred from the iGharv fc Elbazl (|2001l l tem¬ 
plates. 

The IR lumin osity is converted to IR (dust-obscured) 
SFRs using the IKennicuttl (11998H relation. The sum of 
IR- and UV-inferred SFRs are then used to estimate the 
bolometric SFRs. 

5. DUST ATTENUATION OF THE NEBULAR REGIONS 
AND STELLAR CONTINUUM 

SFRs inferred from the UV luminosity and the Ha 
fluxes are shown in the left panel of Figure 0] Without 
dust corrections, the UV and Ha SFRs have a large dis¬ 
crepancy due to the smaller dust absorption cross-section 
at 6564 A relative to that at 1700 A (see the inset panel 
in Figure 0J. The main left panel of Figure 0] shows Ha 
and UV SFRs, both corrected for dust attenuation based 
on the Calzetti attenuation curve and assuming that the 
same color-excess, E(B — V), as derived from the best- 
fit SED model, applies to the stellar continuum and the 
nebular regions. To better quantify the mean trend, mea¬ 
surements have also been performed in bins of UV dust- 
corrected SFR (filled circles in Figure 0|). The bins are 
Alog(SFR[/y )Corr . ecte( i) = 0.3 dex wide. For all the galax¬ 
ies in each bin, regardless of their redsliift, we calculate 
the median flux densities in the observed filters. With 
the new set of median fluxes, the best-fit stellar popu- 
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Figure 4. (Left:) Comparison of SFR(Ha) and SFR(UV), both corrected for dust assuming the Calzetti attenuation curve and assuming 
that the same color excess applies for the stellar and nebular regions. The uncorrected values are shown in the inset panel. Green 
diamonds denote “detected” objects, those whose R s -band magnitudes are more than 1 a brighter than the continuum. Blue arrows 
show the remaining “undetected” objects. Red filled circles are the results of stacked SEDs in bins of dust-corrected SFR(UV), including 
both the detected and undetected objects. Undetected objects bring the stacks lower than the average of the individually detected green 
diamonds. The dashed lines indicate the one-to-one relationships. (Right:) Ratio of SFR(Ha) co7 . T . ecte ^ to SFR(UVj cor .7-ected in bins of 
SFR(UV) cor .rected- The colors denote different recipes for the nebular dust correction. Using the same color excess for nebular extinction 
together with the Calzetti curve (model D; see text) gives the best agreement between Ha and UV SFRs, while assuming higher nebular 
attenuation (i.e., E(B — V) ne b u i ar = 2.27 E(B — V) stellar) with the Calzetti curve results in SFR(Ha) that is significantly larger than 
SFR(UV) (model A). The details of each model are summarized in Table [2] 


lation model is determined through x 2 minimization for 
the mean redshift of the galaxies in each bin. The Her 
and UV luminosities are then derived from the stacked 
SED in the same way as for individual galaxies (see § 0. 
Reported errors are calculated based on the 1 er standard 
error of the mean of the K s magnitudes contributing to 
each bin. 

The appropriate dust corrections to apply to the neb¬ 
ular emission lines is still a subject of debate, as dis¬ 
cussed in Section [I] The attenuation of the nebular lines 
and the UV stellar continuum are not completely de¬ 
coupled as both arise from dust around massive stars. 
On the other hand, whether the color-exc ess is the 
same f or the nebular and stellar regions ( e.g.. lErb et al. 
20064 I Reddy et all 1201 011 or not (e.g.. ICalzetti et ah 
20001: iForster Schreiber et al.ll2009t iGarn et al.ll2010ll ~in 

h igh-redshift galaxies has yet to be fully investigated. 

ICalzetti et al.l (200011 found that E(B — V) ste i = 0.44 x 
E{B — V)neb for a sample of local star-forming galaxies. 
In the absence of direct measurements of the nebular 
color-excess (e.g., via the Balmer decrement), this rela¬ 
tionship can in principle be used to estimate the nebular 
color excess and apply a dust correction to the Ha line; 

e-g-, 

A-Ha — E(B U)nebufar 

= 2.27 K Gal (Ha) E(B - V) stellar 
= 2.27 x 2.52 E(B - V) steUar 
= 5.72 E(B - U)stelfar, (4) 

where /c Gal is the iCardelli et al.l (19891) Galactic extinc¬ 
tion curve, assuming that the nebular regions in the 


high-redshift galaxies abide by such a dust curve. Ajg a 
is the absolute extinction of the Ha emission line, and 
E(B — V) stellar is the SED-inferred stellar color excess. 

We investigate several of the more commonly used 
dust-correction recipes, as described below. The Ha 
SFRs are corrected using five different methods: 

A. the Calzetti attenuation curve for both gas and 
stars, but 2.27x larger color-excess for the nebu¬ 
lar lines than the stellar continuum (as we call it 
“different E(B - V)”), 

B. the Calzetti atten uation curve fo r the s tellar con¬ 
tinuum, and the ICardelli et al.l (1 9891 1 Galactic 
curve for the nebular lines assuming 2.27x larger 
color-excess, 

C. the Calzetti attenuation curve for the stel- 
lar continu um, and an SMC extinction curve 
(Cordon et al.l 1200 3) for the nebular lines assum¬ 
ing 2.27x larger color-excess, 

D. the Calzetti attenuation curve for both gas and 
stars, and the same color-excess (as we call it 
“equal E{B-V)”), and 

E. the Calzetti attenuation curve for th e ste l lar con¬ 
tinuu m, an SMC extinction curve (Gordon et al.l 
l2003ll for the nebular lines assuming the same color- 
excess. 

We choose to compare these different dust-correction sce¬ 
narios for the stacked values instead of the individual 
galaxies as the former represent an average over many 
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Table 2 

Dust Corrected Ha SFRs Using Different Attenuation Recipes, in Bins of SFR(UV) corr - 


SFR(UV) cor . r a 

N b 

Calzetti 

Different E(B — V) d 
(Model A) 

Galactic 

Different E(B — V) e 
(Model B) 

SFR(Ho) cor - r c 

SMC 

Different E{B - V ) f 
(Model C) 

Calzetti 

Equal E(B - V) e 
(Model D) 

Galactic 

Equal E(B - V) h 
(Model E) 

18 ±3 

26 

44 ±3 

38 ± 3 

34 ±2 

30 ±2 

28 ±2 

32 ± 6 

63 

70 ±4 

54 ±3 

46 ±3 

38 ±2 

35 ±2 

65 ± 12 

83 

103 ± 7 

75 ±5 

61 ±4 

49 ±3 

43 ±3 

120 ±22 

56 

200 ± 18 

140 ± 12 

112 ± 10 

88 ±8 

76 ±7 

258 ± 94 

15 

443 ± 88 

295 ± 59 

228 ± 45 

174 ± 35 

146 ± 29 


a The mean and standard deviation of the corrected SFR(UV) in each bin, in units of Mq yr 1 . 
b Number of objects in each bin 

c The SFR(Ha) corr is derived based on the stacked SEDs in bins of dust-corrected SFR(UV), in units of M© yr —1 . The errors are la 
standard error of the mean of the K s magnitudes contributing to each bin, which is converted to the error in SFR. 

d A-Ha — 2.27 K Cstlz (Ha) E(B — V%t e ii a r 5 A# a is the total nebular extinction at 6564 A and K Cstlz (Ha) is the Calzetti reddening at 
6564 A. 

e Aif a = 2.27 K Gal (Fa) E(B — U) s t e iiar 5 ^Gai(^ a ) i s th e Cardelli Galactic reddening at 6564 A. 
f Ah a. = 2.27 KsMc(^ a ) E(B — U) s t e iiar ; K sMc(^ a ) is an SMC reddening at 6564 A. 
s Ah a = ^Caiz (-Ha) E(B — V%t e iiar 5 ^Caiz (-Ha) is the Calzetti reddening at 6564 A. 
h Ah a — K Ga.i(H&) E(B — V’) s teiiar 5 K Gai(^ a ) th e Cardelli Galactic reddening at 6564 A. 

points including both detected and undetected quanti¬ 
ties, each of which may be relatively uncertain. The re¬ 
sults are plotted in the right panel of Figure 00 and are 
reported in Table [5J 

To summarize, these are the values that are used to 
correct the observed Ha SFRs in each model according 
to Equation [U 

a£> = 2.27 K Cala (ifa) E(B - V) = 7.56 E(B - V ), 

A^ } = 2.27 K G jHa) E(B -V) = 5.72 E{B - V), 

A { h1 = 2-27 K SMC (Ha) E{B -V) = 4.54 E(B - V), 

A= «ca dHa) E(B -V) = 3.33 E(B - V), 

A E(B - V) = 2.52 E(B - V). (5) 

In these equations, E(B — V) is the SED-inferred color- 
excess observed for the stellar continuum and A# a is the 
total nebular extinction at 6564 A. 

In model (A) with a larger color excess for the nebular 
regions with the Calzetti curve, Ha-inferred SFRs are 
significantly larger than UV SFRs. Taking the average 
of the SFR(Ha) (model A) to the corrected SFR(UV) in 
the five bins of SFR(UV) indicated in Tables [5] yields a 
value of 2.1 ± 0.2 which shows a ~ 5<r discrepancy from 
unity. 

Using either model (C) or (D), results in 
SFR corrected (TJa) that are in good agreement with 
SFR cor . ree t ec j(C/V). Model (B) reproduce systematically 
higher SFRs but is still consistent with the range 
of SFR cor . rec ted(t7U) in higher bins, and model (E) 
estimations are lower than expected. 

Based on the data provided in Table [2] and the right 
panel of Figurc[4l model (D) provides the best agreement 
between SFRs among the five models. The average of the 
SFR(Ha)-to-SFR(UV) ratios in the five SFR(UV) bins 
of this model shows less than 1 <t discrepancy from unity. 

This suggests that on average the best recipe to correct 
the Ha-inferred quantities of star-forming galaxies at z ~ 

2 is to use E(B - U) ne buiar = E(B - V) ste iiar with the 
Calzetti reddening curve. 



SFR(UV) + SFR(IR) [M 0 /yr] 


Figure 5. The observed Ha line luminosity as a function of 
bolometric SFR. The stars indicate the average values in bins of 
SFR(SED). Galaxies with SFR(SED) < 32 Mq yr —1 are in the 
first bin; those with SFR(SED) between 32 M© yr —1 and 100 
M© yr —1 are in the middle bin; and last bin consists of galaxies 
with SFR(SED) > 100 M© yr -1 . The error in SFR(UV)+SFR(IR) 
is determined by the dispersion error in IR luminosity of the objects 
contributing to each stack. L(IR) dispersion error is estimated 
through bootstrap resampling simulations. The solid blue line is 
the best linear fit. The slope of this fit gives us the conversion factor 
for L(Ha) 0 b s -to-SFR£ 0 £ a z relation. The shaded area is the 68% 
confidence interval of our estimated slope. The orange and dark - 
rcd dashe d lines show t he conversion factors of iKennicutt! (1998) 
and IMadau et ahl (|1998l j. respectively. 

6. BOLOMETRIC STAR-FORMATION RATES 

As discussed above, a primary disadvantage of using 
Ha and UV SFR diagnostics is that we must account 
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for dust extinction, though this is less of a problem for 
1 1 a as it is for the U V. The attenuation curves (e.g., the 
iCalzetti et al.!l2000l curve) that are used to correct the lu¬ 
minosities encode information regarding the dust grain 
size distribution and the geometrical distribution of dust 
with respect to stars (e.g., whether the dust is clumpy or 
uniform, a homogeneous mixture or a foreground screen, 
etc.). On the other hand, the IR luminosity is a di¬ 
rect tracer of dust and thus provides an independent and 
more robust diagnostic of dust attenuation^ Adding the 
IR luminosity, which accounts for the obscured star for¬ 
mation, to the unobscured (UV and Ha) tracers of SFR 
can give a more reliable estimate of the bolometric SFR. 

6.1. The Relationship between Observed Ha Luminosity 
and Total SFR 

To convert the observed Ha luminosity to total _SFR, 
t he m ost commonly used calibration is that of lKennicutti 
(1998;), which is computed using an evolutionary syn¬ 
thesis model, assuming Case B recombination with 
T e =10,000K. In the absence of Balmer decrement mea¬ 
surements, the observed Ha SFR must be corrected for 
dust with some assumption of the attenuation curve and 
color-excess of the ionized gas. In this section, we ex¬ 
amine the Lff a -to-SFR conversion by relying only on the 
observed data, without making any assumptions about 
the dust correction or electron temperature of the ion¬ 
ized gas. With this conversion, one can use the Ha ob¬ 
served luminosity as a proxy for total SFR without as¬ 
suming a dust attenuation curve and its associated un¬ 
certainties. The conversion applies to stellar populations 
with constant star formation at least over 100 Myr. For 
bursty star formation histories over timescales shorter 
than 100 Myr, the Ha and UV+IR luminosities diverge 
as the Ha luminosity traces more instantaneous star for¬ 
mation over timescales of ~ 10 Myr, while IR+UV is not 
sensitive to the change of star formation on timescales 
< 100 Myr. Here, we are deriving the conversion fac¬ 
tor based on the stacked data that represents the “aver¬ 
age” quantities and thus should be used with caution for 
galaxies that may be undergoing bursty star formation. 

Bolometric SFR (SFR(IR)+SFR(UV)) is plotted as a 
function of the Ha luminosity in Figure [5] As both the 
Ha and IR measurements have non-detections, we de¬ 
cided to bin the data with respect to the SED-inferred 
SFR. We chose the bins so that in each bin the SFR(SED) 
is consistent with the stacked SFR(IR)+SFR(UV) within 
the uncertainties, because ultimately we are using 
SFR(IR)+SFR(UV) to calculate the Ha luminosity-to- 
SFR conversion factor. MIPS data are stacked in each 
bin and the total IR luminosity is extracted in the same 
way as for individual galaxies. The UV and Ha lumi¬ 
nosities are derived from the stacked SEDs. 

A weighted least squares regression method is used to 
find the slope of the relation; 

SFR bo , (Mq yr” 1 ) = (7.5 ± 1.3) x 1CT 42 L obs {Ha) (ergs' 1 ). 

( 6 ) 

7 The dust heated by the visible light from older stellar popula¬ 
tions (i.e., the cold component) also contributes to the IR luminos¬ 
ity. However, in case of the star-forming galaxies the contribution 
of the UV radiation of massive stars (the w arm com po nent) is ex¬ 
pected to dominate the total IR luminosity (|KennicuttJll99a) . 


This relation can be used to convert the observed Ha 
luminosity in un its of ergs - 1 to th e bolo metric SFR. The 
factor given in IKennicuttl (1199811 and iKennicutt et al.1 
(I1994I1 . to convert observed/intrinsic luminosity to ob¬ 
served /intrin sic SFR, is 7.9 xlO” 42 M 0 yr -1 erg -1 s. 
iMadau et all (|1998H reported a conversion factor of 6.3 
x 10 42 M 0 yr -1 erg -1 s. These factors are computed 
using evolutionary synthesis models and are subject to 
those models’ uncertainties as well as the initial assump¬ 
tions that went into these models, such as IMF, star- 
formation history, and the stellar evolution and atmo¬ 
sphere m odels . For the same galaxy type and assumed 
IMF, IKennicuttl (1199811 reports ~ 30% variation among 
different calibrations, which mainly reflects the sensitiv¬ 
ity to the SED modeling stellar evolution input data. 
In our analysis, converting L(IR) and L(UV) to SFR 
still rely on the model assumptions, but previous studies 
have shown that in the absence of the Balmer decrement 
measurements, the sum of SFR(IR) and SFR (UV) is 
the most reliable estimate of total SFR (iHirashit a et al.1 
120031 iReddv et all l2012bll . Ou r derived s lope for a 
Salpeter IMF is co nsistent with IKennicuttl (11998( 1 and 
IMadau et all (I1998H within the errors. 

Excluding galaxies with ages < 100 Myr (see § 16.211 
results in a best-fit slope of 7.8 xlO~ 42 M 0 yr -1 erg -1 s 
that is consistent with the slope found above. 

6.2. Validity of the Calzetti Dust Attenuation for Young 

Galaxies 

Figure [5] shows the bolometric SFR vs. median SFR in¬ 
ferred from the best-fit SED model, in bins of SFR(SED). 
Both quantities are representative of the total SFR of the 
galaxies, the only difference being the method by which 
dust is accounted for. The bolometric SFR, as mentioned 
before, is the sum of SFR(IR) and SFR(UV), and is in¬ 
dependent of any dust correction. On the other hand, for 
SFR(SED), the effect of dust is considered by using the 
Calzetti attenuation curve in the SED fitting procedure. 
We would like to determine whether for all types of galax¬ 
ies at z ~ 2 the bolometric SFRs agree with SED SFRs 
that are corrected by the Calzetti attenuation curve. We 
stack the MIPS data of individual detected and unde¬ 
tected galaxies in three bins of SFR(SED) as described 
in ( 12.31 The three bins in log(SFR(SED)) are < 1.5, 1.5 
to 2, and > 2, with respectively 44, 56, and 15 galaxies 
in each bin. SFR(SED) and SFR(UV) are the median 
SFRs of the individual galaxies in each bin. The results 
of stacking is shown with open blue circles in Figure [6l 
In the last two bins, the SED SFR is overestimated com¬ 
pared to the bolometric SFR. Stacking the MIPS images 
for only objects that are inferred to be older than 100 Myr 
produces different results. The number of objects in bins 
with galaxies older than 100 Myr reduces to 43, 50, and 8, 
respectively from the lowest SFR(SED) bin to the high¬ 
est. 

Once we remove galaxies with ages < 100 Myr, 
SFR(SED) agrees well with SFR(UV)+SFR(IR). An 
age-dependent L(UV)-to-SFR calibration (as we use 
here) does not decrease the difference between 
SFR(SED) and the bolometric SFR for stacks of all 
galaxies. It seems younger galaxies are at the same time 
less dusty (smaller bolometric SFR) but redder (larger 
SFR(SED)) than their older counterparts. These re- 
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Figure 6. Comparison of the bolometric SFR and SED-inferred 
SFR. The values are divided into three bins of SFR(SED) includ¬ 
ing individual detections and individually undetected galaxies at 
24 /zm; Open blue circles include all objects and filled red circles 
include only those galaxies older than 100 Myr. Diamonds are indi¬ 
vidual detected objects and downward arrows represent 3cr upper 
limits in SFR(IR)+SFR(UV) for the MIPS undetected galaxies. 
Young galaxies (ages < 100 Myr) are colored blue and galaxies 
with ages > 100 Myr are orange. SFR(IR) error bars are estimated 
through bootstrap simulations and show the dispersion in the SFRs 
contributing to each stack. The last bin consists of galaxies with 
SFR(SED) > 100 Mq yr — by removing young galaxies (ages < 
100 Myr) about half of the sample in this bin is dismissed and we 
see the change between the two measurements (including and ex¬ 
cluding the young ones). The dotted line indicates a one-to-one 
correspondence. 

suits are con s istent with those of I Reddv e t, all (12006H1 . 
IReddv et alJ |2010), and IWuvts et al.l (|2012D . who find 
that galaxies with young stellar population ages are less 
dusty for a given UV slope than older galaxies. An 
SMC-like curve may be more appropriate for these galax¬ 
ies. One possibility is that young galaxies have la rger 
dust covering fractions than older galaxies (IReddv et all 
120101 1. A large dust covering fraction makes the UV slope 
redder for a given amount of dust attenuation. This 
study confirms that at z ~ 2 the Calzetti attenuation 
curve is applicable to star-forming galaxies older than 
100 Myr, but a steeper attenuation curve may be neces¬ 
sary for the younger galaxies. 

7. SUMMARY 

We have studied the multi-wavelength properties of 
a sample of ~ 200 galaxies with spectroscopically con¬ 
firmed redshifts at 2.08 < 2 < 2.51, to study the validity 
of commonly used dust correction recipes and to compare 
SFRs inferred from the Ha line luminosity, the UV con¬ 
tinuum luminosity, and Spitzer/MIPS 24 /jm measure¬ 
ments that are converted to the total IR luminosity. In 
this study, we benefit from using the broad-band photom¬ 
etry excess techniques to determine the Ha fluxes and, 


hence, conduct a large sample of 2 ~ 2 Ha measurements 
that is immune to uncertainties in the spectroscopic slit- 
loss corrections. The galaxies’ properties are determined 
from stellar population model fitting to the rest-frame 
UV through near-infrared data. The main conclusions 
are as follows. 

• By investigating several recipes for dust correct¬ 
ing the nebular emission lines, we find that assum¬ 
ing the same color excess of the ionized gas and 
the stellar continuum (i.e., E{B — V) ne b = E(B — 
R)stei), and assuming that the Calzetti attenuation 
curve applies to both, results in the best agree¬ 
ment between SFR(Ha) and SFR(UV). If we as¬ 
sume the Calzetti curve to both the stellar and neb¬ 
ular regions but use E(B — V) ne b u iar=2-‘27 E(B — 
V)stellar, the corrected SFR(Ha) measurements 
are inconsistent with the corrected SFR(UV)s at 
~ 5cr level, averaged on the five bins of SFR(UV) 
(see the right panel of Figure HJ). 

• Using the available Spitzer/MIPS data for ~ 100 
galaxies in our sample, we derive an observed L Ha~ 
to-SFR totai conversion factor of (7.5 ± 1.3) x 10 -42 
Mq yr _1 erg -1 s. This calibration is independent 
of any assumptions on the dust correction and can 
be used to convert the observed (extincted) Ha lu¬ 
minosity to a bolometric SFR when no dust at¬ 
tenuation measurements for the Ha luminosity is 
available. 

• By comparing the stacks of SFR(UV)+SFR(IR) 
with the SED-inferred SFRs that are corrected 
for dust by the locally-derived Calzetti curve we 
conclude that a steeper attenuation curve (such 
as an SMC curve) may be necessary for galax- 
ies younger than 1 00 Myr, as previous st udies 
(|Reddv et alJl2006bL 1 201 di Wuvts et al. 1120121) have 
found. We find that applying the Calzetti curve to 
the stacks of all galaxies, including the young ones, 
results in SFR(SED)s that are inconsistent with 
the SFR(UV)+SFR(IR) at ~ 2cr significance (the 
SFR(SED) in three mass bins of Figure [6] are incon¬ 
sistent with SFR(UV)+SFR(IR) by 1, 2, and 3cr, 
from the lowest to highest mass bins respectively). 
The young galaxies have redder UV slope and at 
the same time lower bolometric SFR compared to 
their older counterparts. As a result, applying the 
Calzetti curve to the young galaxies overestimates 
the SED-inferred SFR when compared with their 
bolometric SFR. The Calzetti attenuation curve 
shows a good overall agreement for galaxies older 
than 100 Myr. 

A more detailed investigation will include direct tracers 
of nebular dust extinction (i.e., the Balmer decrement). 
In the future studies, we plan to investigate this aspect 
with the MOSFIRE Deep Evolution Field (MOSDEF) 
survey, which uses the near-IR multi-object spectrograph 
MOSFIRE on the Keck I telescope to obtain spectro¬ 
scopic measurements of the nebul ar emissio n lines for a 
sample of ~ 1500 galaxies dKriek et al. 1120141) . 
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